UCLA/92/TEP/41 



Explicit Calculation of the Renormalized 
Singlet Axial Anomaly 



Michiel Bos* 

Department of Physics 
University of California 
Los Angeles, CA 90024 

A careful and complete discussion is given of the renormalization of the singlet ax- 
ial anomaly equation in a vector-like nonabelian gauge theory such as QCD regularized 
by ordinary dimensional regularization. Pseudotensorial structures are treated with the 
't Hooft-Veltman prescription. A general framework for calculations is developed, and 
subsequently verified by explicit computations through two loops. This is followed by a 
discussion of the matrix elements obtained. 



11/92 

* e-mail: bos@uclahep.jnet, or bos@uclaliep.pliysics.ucla.edu. 



1. Introduction 

The axial anomaly is by now a standard feature of gauge field theory. 

In spite of this the study of its renormalization has still not reached a state of 
stable equilibrium. The literature on the question, which began more than twenty 
years ago with the classical work of Adler and Bardeen has grown both by 
the emergence of new issues and by the reconsideration of old ones, a process that 
continues to the present day. Particularly notable revivals of interest have been 
spurred by the paradox of the currents in supersymmetric theories (see e.g. 
for an early formulation of the problem) and, more recently, by new developments 
concerning the spin of the proton (reviewed e.g. in 0). 

The present paper intends to contribute to the topic an explicit perturbative 
computation at the two-loop level for a vector-like nonabelian gauge theory such 
as QCD. More precisely, I will calculate the two-gluon matrix elements of all the 
operators in the singlet anomaly equation to second nontrivial order in dimension- 
ally regularized perturbation theory, and apply the results to an analysis of the 
equation itself. The choice of ordinary dimensional regularization (as against the 
dimensional reduction of much of the supersymmetric literature, e.g. 0) is based 
on its prominence in perturbative studies of nonabelian gauge theories, as well as 
its relatively obvious consistency. The prescription of 't Hooft and Veltman is 
adopted for the pseudotensorial objects 75 and e^^P'^ . This type of approach is of 
course not novel in itself; the present investigation is similar in spirit to several 



earlier ones, notably ^ and . I have gone beyond this earlier work by perform- 
ing calculations that are both more comprehensive and more complete; I have also 
gone to greater lengths in connecting the actual computations to their conceptual 
underpinnings. 

The first two sections following this introduction are devoted to a general 
derivation of the renormalized anomaly equation. First I review the operators of 
interest, normalizing them with reference to their scaling properties to the extent 
allowed by universality; renormalization group invariance of a renormalized anom- 
aly equation then leads to strong a posteriori constraints on its coefficients. This 
line of argument, which originated in early work on the intertwining of the anomaly 



and Callan-Symanzik equations [11,12|, has by now become rather commonplace 
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in the literature, where one finds it expounded with varying degrees of precision. 
Next I adopt a constructive approach, regularizing the theory via dimensional reg- 
ularization, deriving the anomaly equation and renormalizing it, and matching the 
result with the scheme-independent template obtained before. 

Section 4, which is the core of the paper, contains most of the actual calcula- 
tions performed. Apart from the very validity of the renormalized anomaly equation 
they verify the finite renormalizations required in the formalism of dimensional reg- 
ularization. In view of the role of renormalization group arguments in organizing the 
problem I have employed the background field method, within which the renormal- 
ization of the nonabelian theory is most analogous to that of the well-understood 
abelian case. I pay particular attention to the normalization of the gluon tensor, 
which has been taken for granted in previous treatments along similar lines. 

The final section is devoted to comments on the matrix elements obtained, 
mostly from points of view suggested by the abelian Adler-Bardeen theorem 
The discussion eventually touches on an ongoing discussion about a possible renor- 



malization of the vacuum angle in QCD |T3,14,15 



The literature on anomalies is so extensive as to preclude exhaustive citation. 
My list of references is a selection of the papers best known to me, and I apologize 
to authors of omitted work. 

The system considered consists of a generically nonabelian gauge field A with 
coupling constant (7^, and fermions '0 of mass m forming rin copies of an irreducible 
representation R of the gauge group. The group invariants cn and Tji are defined 
in the usual manner, tr T"T^ = TRdab, T^T" = crI, where the T" are the group 
generators, [T",T^] = ifabcT'^] in particular, for the adjoint representation CASab = 
facdfbcd- I denote the dimensionality of the regulating spacetime of dimensional 
regularization by d and write d = 4 — e. 



2. Operator Normalization 

This section reviews the general structure of the anomaly equation. The theory 
is assumed to be regularized in some gauge-invariant way, but the details of the 
regularization will be irrelevant. The point is to describe a sensible normalization of 
relevant operators by concentrating on scaling behavior and universality arguments. 
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The following auxiliary consideration will be required. Suppose one is given 
a renormalized (finite) operator O; suppose moreover that its renormalization was 
multiplicative and its anomalous scaling therefore of the form 

fi^O = -^o{9')0 (2.1) 

where is the coupling constant and is a renormalization mass. (I will write all 
expressions in such a way that they are naturally interpreted in dimensional regu- 
larization with minimal or quasiminimal subtraction by taking fi to be the familiar 
scale factor of the extra dimensions, but no firm commitment to that procedure is 
necessary.) A finite multiplicative renormalization of C is a replacement of O by 
O' = f{g'^)Oj fig"^) being a power series in g"^ with /(O) = 1; the latter restriction 
excludes redefinitions of the classical operator corresponding to O that have noth- 
ing to do with the ambiguities of renormalized perturbation theory. The anomalous 
dimension of O' differs in general from that of O and is given by 

lO'{9')^lo{9')-(3{9')j^. (2.2) 

The power series for the beta-function starts at order g"^ and that of 70 typically 
at order g'^, but if 70 has no one- loop term, 70 = cg^ + the quotient 70 //3 
is a power series in g^ and a finite renormalization / may be found such that the 
renormalized operator O' has vanishing anomalous dimension, 

/te^)=exp[/°'d.||)]. (2.3) 

On the other hand, the part of the anomalous dimension that arises at one loop 
is unaffected by finite renormalizations, hence irremovable. Similar reasoning ap- 
plies to additive renormalizations; terms of order g"^ and beyond in the power series 
expansion of an off-diagonal anomalous dimension are removable by a finite renor- 
malization, while the one-loop, g^ term is universal. 

Anomalous Ward identities for the singlet axial current relate insertions in 
Green functions of three operators: the divergence of J5 , the singlet axial current 
itself; j5, the pseudoscalar density; and the anomalous term FF. In addition they 
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involve contact terms corresponding to axial transformations of the fermion fields. 
To lowest nontrivial order one finds the familiar expression 

d,J^ = 2mj5 + ^^^^e'^^^'^F;,^;, + c.t., (2.4) 
"c.t." being the contact terms. Upon insertion in a Green function of fermion fields 
'ip{xi), ...,ip{xp), 'ip{yi), ...,'ip{yq) and gauge fields A{zi), A{zr), ''^c.t." takes the 
form 

p 

- ^S{w - Xi) < ijj{xi)...'y5^jj{xi)...^jj{xp)ijj{yi)...ij{yq)A{zi)...A{zr) > 

Q 

- ^ 5(W - l/j) < ij{xi)...ij{Xp)'lJj{yi)...'lJj{yj)l5--4{yq)Mzi)---MZr) > 
J = l 

(2.5) 

involving the axial transformations of the fermions. 

Classically the singlet axial current equals 'ip^^'jsip'i its quantum version is 
renormalized multiplicatively. The one-loop diagram contributing to that (Figure 
1) is identical to the diagram for the vector current 'ip^^^'tp but for an extra 75 
at the vertex. The inclusion of this 75 does not spoil the familiar cancellation of 
infinities between the one-loop correction and the wave function renormalization of 
the constituent fermion fields, so J5 , too, is finite at one loop and has no one- loop 
anomalous dimension. By the argument reviewed above one infers the existence of 
a (unique) normalization that makes the scaling of the axial current canonical to all 
orders; call this current [jsjc- We may and will take the renormalized divergence 
of the axial current to be the divergence of the renormalized axial current. 

Similar reasoning applies to the pseudoscalar density J5, classically equal to 
^V'TsV'- Its multiplicative renormalization involves at one loop a diagram as in 
Figure 1. Now an analogy obtains with the one- loop renormalization of t/St/^, the 
renormalization constant of which is the inverse of that of the mass. The combi- 
nation mj5 is thus found to scale canonically at one loop and may be taken, via a 
finite renormalization of if necessary, to have that property to all orders. I denote 
the resulting object by m[j^]c- 

The composite operator FF = e'^'^'"^ F^^^F^^ is the derivative of the gauge- 
variant Chern-Simons current Ae^'''P'' [A^dpAl + \fahcAlA^pAl). The renormaliza- 
tion of this current is not exhausted by a multiplicative factor, but involves mixing 
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with the axial fermion current as well. (By contrast, the axial current itself is 
protected by its gauge invariance from non-diagonal renormalization involving the 
gauge-variant Chern-Simons current.) Correspondingly FF mixes with d^j^. 

Figures 2 and 3 contain one-loop diagrams relevant for the diagonal and off- 
diagonal renormalization of FF, respectively. The analysis of Figure 2 is simplified 
if background field Feynman rules are used, and I will do so from now on. The 
logarithmically divergent parts of the diagrams in Figure 2 then cancel among each 
other, so the diagonal, multiplicative scaling behavior of FF is derived entirely 
(at least to this order) from the wave function renormalizations of its constituent 
gluons. In the background field formalism the renormalization of the gluon field is 
by construction exactly the inverse of that of the coupling constant, so the product 
g^FF is not renormalized and free of anomalous dimension at one loop, again as 
far as diagonal renormalization is concerned. The option of a finite renormalization 
of FF may then be invoked to extend this property to all orders. A renormalized 
version of FF will be denoted by [FF]c if its product with g'^ has canonical diagonal 
scaling. 

This leaves the issue of additive renormalization and off-diagonal scaling, to be 
parametrized by 

i^^i9"[FF]c) = -iF,{g') g'd,\j^]c. (2.6) 

The extra g^ on the right hand side balances the one on the left in order that the 
counting of powers of g^ in the anomalous dimension be identical to that of 
the number of loops at which they arise. Evaluation of the diagram in Figure 3 
shows 7f5 to be nonvanishing at one loop; the logarithmically divergent part of 
that diagram is of the form 

in momentum space, A being a cutoff and jj. a renormalization scale, and p and 
p' the momenta of the fermions. This implies a one-loop value ^pb = 3cij^^/27r^. 
Higher orders may be changed by additive renormalizations of [FF]ci {[FF]c)' = 
[FF]c + f{g'^) dix[i^]c with /(O) = 0, but there is no particularly natural way 
of fixing this freedom and I will not commit myself with regard to it, leaving the 
ambiguity of [FF]c implicit in the notation. 
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(There exists another approach to the normahzation of g^FF, ultimately equiv- 
alent to the present one, in which the Chern-Simons current is related by gauge or 
BRS descent to a naturally finite operator, the normalization of which fixes the 
diagonal normalization of the renormalized g^FF; see [1T6| , [T7| , |T8| . ) 

Now assume that the anomaly equation survives renormalization, in the sense 
that there exists in the complete theory a relation of linear dependence involving the 



renormalized operators defined above and reducing to ( |2.4| ) at the lowest nontrivial 
order in perturbation theory. Write this relation in the form 



Ci{g') d,[j^]c = C2(/) 2m[j,]c + Cs{g') ^/[FF]c + c.t.. (2.8) 

where the Ci are power series in (7^, with Ci(0) = 1 to accommodate (|2.4|), and 
the overall factor has been fixed by normalizing the contact terms to their standard 



value, (|2.5|) . Application of the renormalization group operator to both sides 
yields 

PC[d,[j^]c = f3C!,2m[js]c + f3C^^^g'[FF]c-Cs'^g'^F6d,[j^^^ (2.9) 

The contact terms have disappeared because they are finite as they stand, see ( |2.5|) , 
hence without anomalous dimension. Without them the remaining operators are 
independent and consistency requires 

Q = C5 = 0, Ci = - "f^/';" C3. (2.10) 

i.e. 

C, = l-'^ r\.nis(£>, C. = C. = 1. (2.11) 

16n^ Jo /3{x) 

The single remaining non-universal feature (at least by the criteria espoused in the 
above) may be hidden in notation by introducing a finite renormalization of the 
canonical axial current, 

[JnAB^C^{g'M]c, (2.12) 
in terms of which the renormalized anomaly equation reads 

d,[j^UB = 2m[j,]c + '^g'[FF]c + c.t. (2.13) 
In keeping with the (mostly supersymmetric) literature [JsIab may be called an 
Adler-Bardeen current, with indefinite article to signal the normalization ambi- 
guity which it inherits, via Ci, from the mixing dimension 7^5- Along comes a 
multiplicative anomalous dimension starting at order g'^ @,^|- 
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3. Dimensional Regular izat ion 

The next step is to derive the anomaly in the framework of dimensional regular- 
ization. This discussion is complementary to the previous one in concentrating on 
aspects that are peculiar to a specific regularization scheme. For convenience I will 
mostly adhere to the conventions and notation of Collins' textbook on renormal- 



ization [20], where the Ward identity for the (non-anomalous) non-singlet current 
is derived. 

To define the axial current in the regulated theory requires a definition for 75, 
which is notoriously problematic in dimensional regularization. In the 't Hooft- 
Veltman prescription 75 is taken as the product of the first four 7-matrices 



at the cost of anticommutativity with all 7's in d dimensions. This requires a 
distinction between the first four dimensions and the remaining (d — 4), to which 
end one introduces projectors denoted by a superscripted bar and hat, respectively; 
a vector is written as d^ -|- a^, with d^ = if n = 0, 1, 2, 3 and otherwise, and 
the other way around for a^. One may then proceed to write down identities for 
tensorial objects with contracted indices. The totally antisymmetric tensor e^i/po- is 
introduced as a bar-carrying object in the above sense, i.e. e^^p^ equals its usual 
value when each index is 0,1,2, or 3, and otherwise; for brevity of notation the 
bar will not be indicated explicitly. Finally, 75 is defined by 

75 = ^e,.,.rYYr- (3.1) 
Its most notable feature is its failure to anticommute with the hatted 7's: 

{75,7/^} = but [7^,71 = 0. (3.2) 

Trace identities may be established as usual: the trace of 75 times a product of n 
7-matrices vanishes if n is odd or less than four, while (taking tr 1 = 4) 

tr ^^^^'YY1'' = ^ie^'"^"- (3.3) 

The bare singlet axial current and pseudoscalar density in the dimensionally 
regulated theory are now defined by 



35 = ^o7''75V'o 
7 



(3.4) 



is = ^^075-00- (3.5) 

ipo = {Z2)^ip is the bare fermion wave function. The bar on in the axial current, 
which does not affect the classical limit, is necessary for hermiticity. 

The would-be conservation law for the axial current is the unrenormalized 
anomaly equation, 

d^j^ = 2moi5 + Oan + 0,t (3.6) 
where Oan and Oct are given by 

a, = iV'o{¥,75}^o (3.7) 

and 

Oct = itpo'j5ii^ - "^o)'0o - #o(^^+ "^0)75^0- (3.8) 

The operator Oan represents the potential anomaly. Only the hatted components 
of the covariant derivation occur; Oan is "evanescent" and would vanish if 
no divergences occurred in the theory in the limit d ^ 4. Oct is proportional to 
the equation of motion and vanishes on shell, but its insertion in covariant Green 
functions, where derivatives and time-ordering are taken to commute, produces 
the contact terms associated with axial transformations of the fermion as given 
explicitly by (|2^ ). 

The next step is to renormalize the equation, for which we adopt minimal 
subtraction. The linearity of that procedure implies that the subtracted equation 
is obtained by replacing each operator by its minimally subtracted version (and 
mo by the minimally renormalized mass m; we assume the coupling and wave 
function renormalizations such as Z2 to be defined by minimal subtraction as well). 
The minimally subtracted current [is ]ms and pseudoscalar density [j'sIms are both 
multiplicatively related to the respective unrenormalized operators. The contact 
terms are manifestly finite, [Oct] MS = Oct- These relations and the anomaly equa- 
tion implicitly determine how the minimally subtracted anomaly operator [Oan] MS 
is to be defined in terms of unrenormalized operators and renormalization factors. 

The crucial point is now that in the limit d 4 finite contributions involv- 
ing the subtracted anomaly term can only arise in diagrams where the evanescent 
vertex on which that term is built multiplies an internal divergence. As a result, 
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[C^aulivis survives in the limit d—>-4 as a linear combination of local operators with 
non-evanescent basic vertices, the net effect being a reduction in the number of 
independent operators. The coefficients of this linear combination are (at least in 
perturbation theory) power series in without constant terms, because [Oan]MS 
vanishes at the tree level. Other examples of this phenomenon have been worked 
out in contexts rather different from the present one P^ , ^ . 

The set of non-evanescent operators in terms of which [Oan]MS is expressed may 
be restricted by dimensional and symmetry considerations. That leaves ^^[jsjMS, 
"^[JsIms, Oct and [FF]mSj the latter being the minimally subtracted version of 
the composite operator FF = e^'^'P" {F^^)o{F^^)o. ((F^^)o = (^3)^ (i^;?^)o is the 
bare field strength, simply proportional to the renormalized one in the background 
field method.) Od may be replaced by other operators by use of the minimally 
subtracted anomaly equation and the remainder may be rearranged as 

PanUs = h{g'')d^,[3^]MS + /2(^7')m[j5]MS + U9^)[FF]ms (3.9) 

where the fi are all finite and of order g"^ at least. Substitution in the minimally 
subtracted anomaly equation yields 

(1 - h{,g''))d^[j'^]MS = (1 + i/2(^7'))2m[i5]MS + h{g'')[FF]Ms + ^ct- (3.10) 

The existence of the anomaly thus amounts to /a 7^ 0. 

The factors 1 — /i and 1 + ^/2 may be removed by a finite renormalization, 
leading to oversubtracted (with respect to minimal subtraction, that is) operators: 

[f,]os^{l-h{9^)M]MS (3.11) 

[k]os = {l + \f2{9^))[k]MS. (3.12) 

An analogous oversubtr action may be performed on FF; fs, if nonzero, has the 
form f^^^ g^ + 0{g'^) with f^^^ 7^ 0, so f^{g'^) / { f^^ g"^) is a power series starting with 
1 and I renormalize 

[FF]os^^j^AFF]MS (3.13) 
leading to a rewritten renormalized anomaly equation 

dM]os = 2m[j5]o5 + fi'^g'[FF]os + O^t- (3.14) 

9 



This brings the construction to the point where contact is made with the dis- 
cussion of the previous section. This time the renormahzed anomaly equation has 
been derived, rather than assumed. Comparison with ( |2.13|) identifies [j^]os as an 
Adler-Bardeen current and mfjsjos as the canonicaUy scahng pseudoscalar density 
"^[i5]c; 9'^[FF]Qg is a product of the type g'^[FF]c with canonical diagonal scaling 
behavior. The value of the off-diagonal anomalous dimension of g'^[FF]Qg beyond 
one loop as well as the related anomalous dimension of [j^]os emerge as specifics 
of the regularization. 



4. Calculations 

In this section I describe the calculation of quantities relevant to the next-to- 
leading order two-gluon matrix element of the anomaly equation. 

The kinematical setting is as follows. We will compute one-particle irreducible 
functions of two gauge fields with a single operator insertion, that operator being 
one of those that occur in the anomaly equation. The labeling of the external 
gauge fields will be as in Figure 4, with (a, b) and (/x, group and Lorentz indices, 
respectively, and p and —p' the incoming gauge field momenta. For the operators 
of interest O the two-gluon function has the form 

e^^^paP^p'^ Sab f H (4.1) 

where To, to be referred to as the reduced matrix element, is a function of the scalars 
and p ■ p' as well as of the fermion mass m. Where an expression for Fq for 
all p and p' is not practicable or illuminating I will specialize to the infrared regime 
p' ^p, accessible from the complete function after differentiation with respect to p, 
say. I will then also eliminate either p^ or ; in each case the renormahzed reduced 
matrix element, classically dimensionless, will depend on the remaining parameter 
at most via renormalization group logarithms. All computations are performed 
using the background field Feynman rules, a convenient table (and motivation) of 
which may be found in |]2^ . 



The following ingredients are required for an evaluation of the two-gluon matrix 
element of the anomaly equation ( 3.14|) through next-to- leading order. The matrix 
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elements of ^^[jsjMS and mfjsjMS) which have their lowest nonzero contributions 
at one loop, should be calculated at one and two loops. The matrix element of 
[FF]ms starts at the tree level and we need the one- loop contribution in addition 
to that. The functions fi{g^) and f2{g^) contribute to the anomaly equation from 
the next-to-leading level onward and their leading (i.e. g"^) terms suffice. Of the 
function fsig'^), on the other hand, both the g'^ and the g'^ terms are required. The 
two-gluon matrix element of Oct is zero. 

We begin with the reduced matrix elements i'df.lJi^lMs ^rn[j5]Ms- '^^^ lowest 
order contributions stem from the familiar triangle graph. Figure 4, and are easily 
computed for general p and p' . The results have been well-known since the early 
investigations of the anomaly, 

^d^[j^]MsiP^P''^^)\one loop = _ ^ ^^^^ Q + W^/qO (p, / ; w)^ (4.2) 
^m[j5]MsiP^P''^'^)\one loop = _ ^ lQo{p,p' \m) (4.3) 



where 



1 /.l-s 



^0 



ho{p,p'\rn)= I ds dt ^— — — — — -. (4.4) 



p^t{l — t) +p''^s{l — s) — 2p -p'st — m? 



If the fermion mass vanishes, 

27r2 



^dMUsi.P^P'^^^^)\one loop = ^^J"^ , f ^[j-^]^^ (p, p'; m = 0) = 0; (4.5) 



by contrast, if the fermion mass is finite but p'=p,p — >0 

^d^[j^]MsiP=P' ^P"^ = ^'^'^)\one loop = 
^ m[js,]MS\P~P jP —^i'^)\one loop — 

We proceed to the contributions of order g'^. The relation between minimally 
subtracted and unrenormalized reduced matrix elements is 

^[0]ms = ^oZA^OMre (4.7) 

with Zq the renormalization constant of O as determined by minimal subtraction, 
O being any of the operators of interest; Z\ is the wave function renormalization 
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of the gauge field. For arbitrary p, p' and m the expressions for the reduced matrix 
elements are rather unwieldy and I concentrate on the aforementioned limits. 

At the two-loop level i'df.j^ ,bare s-iid Tj^^bare iuvolve the eleven diagrams of 
Figure 5. The values of these diagrams with (1) = d^jt^, P=p', P^T^O, m = 0, (2) 
O = d^j^, p = p' , ^ 0, m 7^ 0, and (3) O = js, p = p' , ^ 0, m 7^ 0, are listed 
in Tables 1, 2 and 3, respectively. They have been computed using the background 
field Feynman rules. As we only determine terms through order g'^ of the left hand 
side of ( |4.7|) , the bare parameters Qq and mo have been replaced in the tables by 
the lowest order of their renormalized expansions, g'^fi'^~'^ and m. A common factor 
g'^riRTu/ (47r)^ has been omitted; Lm stands for In + 7 and Lp for In + 7. 
The gauge dependence of the gluon Feynman rules becomes an issue as one finds 
terms up to quadratic order in the standard covariant gauge parameter a. 

Consider first d^j^ for massless fermions. The two-loop diagrams of Table 1 
add up to 

^d^j^,bareiP = P'^P'^7^0'^^ = 0)\two loops ~ 

g^'riRTR f 3 1 (4-8) 

- g^4 (cfl + c^(l + -(l-a)--(l-a) )j. 

The similarity sign is meant to express the replacement of bare parameters by 
their lowest order renormalized ones. The absence of a term proportional to cr a 
corresponds to the gauge invariance of the result for QED and obtains equally from 
the other two tables. It is not hard to show by diagrammatic reasoning that such 
terms must cancel independently of the values of p' or m. It may also be noted that 
the propagator-induced corrections to the triangle, diagrams I through III, vanish 
at a = 0, a consequence of the vanishing of the one-loop self-energy of massless 
fermions in the Landau gauge. Proceeding to the counterterms, the bare one-loop 
matrix element is _ £ii^^^ii(^p2-j(d-4)/2^ ( |4.5|) with the regulating dimensionality d 
taken into account; but the bare coupling in this expression is converted to the finite 



renormalized one by multiplication with Za as in (|4.7| ), because the renormalization 
of the wave function is precisely the inverse of that of the coupling in the background 
field method. As for the renormalization factor Zq^jI^'-, it is equal to 1 at one loop 
by the calculation sketched in section 2. The renormalized reduced matrix element 
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becomes, through two loops, 



^d,[jt:]MsiP=P ^0;m = 0) = 

If the fermions are taken to be massive the two-loop unrenormalized matrix 
element of S^Js is zero in the limit p = p' , ^0, as is found by adding the entries 
of Table 2; as a result the vanishing of the renormalized matrix element persists to 
order g"^, 

^dAj^Us(P=P'^P'^^-^^7^^)=0- (4-10) 

This extends to the nonabelian case the corresponding property of the divergence 
of the axial current in massive QED, proven to all orders by Adler and Bardeen 
we shall return to this presently. 

Turning to the pseudoscalar density j5, which occurs in the anomaly equation 
multiplied by m and is therefore relevant only in the case of massive fermions, the 
contributions to the limit p=p', p^—^0 tabulated in Table 3 add up to 

'i^j5,bare{P = P',p'^ ^O;'m^0)\two loops ~ 

g^'uRTR fl /I ^3 1 (4-11) 

Note that though the infrared behavior of individual diagrams produces terms in- 
volving logarithms of p'^ , such terms cancel in the sum. The bare one- loop result is 
here 

f,5,ba.e(p=p',p'-0;m^0)|„„e loop = ^^^(^o)'-'. (4.12) 

The counterterms implicit in gg cancel against those generated by Za, as before. A 

similar cancellation arises between — and Z^- ; we saw in section 2 that the one- 
mo -'5 ' 

loop renormalization of is opposite to that of m, mojs = "^[jsIms- That leaves 
only the scale factor, which by use of the familiar minimal subtraction for the mass 



mo = (1 - + . . .) m (4.13) 
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is found to contribute 3g^ciinjiTji/327r'^m to the reduced matrix element. We ob- 
tain, to two loops. 



(4.14) 



Thus f2m[j5]Ais for massive fermions is similar, but not completely identical, to 
^m^sIms massless ones. 

The results above may also be arrived at by subtracting sub divergences diagram 
by diagram. This allows for the special case of the Feynman gauge comparison with 
some of the ordinary perturbation theory results of [|10|, which are presented in 
this form; I have found agreement wherever I checked. Upon further restriction to 
massless fermions one also finds agreement with the older results of . Restriction 
both to massless fermions and to the group structure cr produces a value given in 



The next task is to find the functions fi{g^) and /2(fi'^) at order g'^. This may 
be done by computing the one-loop matrix element of the anomalous operator Oan 
with a fermion and a conjugate fermion field; the tree graph vanishes in the limit 
(i — s> 4. The relevant one-loop diagrams are shown in Figure 6. After inclusion of 
the wave function renormalization of the constituent fermions of Oan one finds the 
finite result, in a general gauge a, 

^^(F -?)75 ^"^75- (4.15) 

This is rewritten in the form of the right hand side of ( p.9|) by comparison with the 
basic vertices of d^j^ and js, equal to i(J — ^)75 and Z75, respectively, and observing 
that the i/j-i/j matrix element of FF starts only at order g^. It follows that 

= + A(9^) = -^ + --. (4.16) 

Obviously this is basically an abelian effect, hence independent of a. 

The above results may now be combined into two-loop expressions for the 
reduced matrix elements of the oversubtracted fermionic operators d^[j^]Qg = (1 — 

14 



fi)d,j.\j^]MS and [j5]os = (1 + |/2)[J5]ms in a general covariant gauge a: 

ra,[,-]^^(p=P,P ^0;m = 0) = 

g^CAURTRf 3 X In ^ 



(4.17) 



fQ^y.,{p=p',p^^O-m^O) = 0. (4.18) 



^m[j,]osiP=P ^0;m^0) = 



47r 



2 



q'^CAnpTB ( 3 , x 1 / 



(4.19) 



The finite renormalization has had the effect of ehminating terms involving the 
quadratic Casimir of the fermion representation from the g'^ contributions. 

We now turn to the matrix element of \FF\ms- The tree level value, determined 
by the basic vertex of FF, is V ^pp^^\^tree = —8. The next order is obtained from 
a formula of the form (|4.7| ). The unsubtracted one-loop diagrams are shown in 
Figure 2; in accordance with the considerations of section 2, their sum is finite if 
is replaced by g'^\i^~'^ . The renormalization constant Zpp must therefore be chosen 
to compensate for the wave function renormalizations of the gluon fields in FF ^ 
Zpp = {Za)~^- With this choice, i'[FF]Ms loop is found to be, for arbitrary 

1 r r^ — s 1c/ 

+ -{l-a) 1+ ds dt ~ ~ {{1 - 2s)p^ + {1 - 2t)p'^ + 2{s + t)p ■ p') 
4 Jo Jo ^ 

(4.20) 

with Iqq as in (^j^) and 

A = p^t{l -t)+ p'^s{l -s)-2p- pst. (4.21) 
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This somewhat unsavory expression has an easily read off hmit as p' ~^p; including 
tree and one-loop effects, 

^[FFUs(P'=P) = -8(1 + + - «) - - (4-22) 

Note that to this order the fermions play no role in '^[pp^j^jc;- 

The one-loop matrix element of FF was computed using a different interme- 



diate regularization in |^^. There the term proportional to (1 — a) is presented 
as 

-£if!(l - - l(p - +p'^5^)/oo] (4.23) 

where the partial derivatives are defined by treating /oo(p, p'', 0) as a function of the 
invariants p^, p'"^ and p ■ p'. This appears to yield ^(1 — a) instead of |(1 — a) in 
( [1.22| ) []TB|, but Iqo is singular as p'— 

dp'^ ^ dp''^' {p' — pY 

and this eliminates the apparent discrepancy (cf. also jT^). 



With an eye on ( ^.141 ) we turn now to the determination of the Taylor 
coefficient of g'^ in fs{g'^). As with /i and /2, the g"^ term in /a may be extracted 
from a one- loop diagram, in this case the triangle with two gluons (Figure 4) . With 
Oan at the upper vertex this graph yields, in the limit (i^4 (to eliminate evanescent 
contributions) Q 

ro,„(p,p';m)|o„e loop = o^2^^ - (4-25) 
This is compatible with (|3.9| ) provided 

/3(/) = ^^ + ..., (4.26) 



as the first two terms on the right hand side of (|3.9| ) contribute to '^Oan oiily 
the g'^ level, is thus determined to be njiTu/ 16%"^, and one recognizes in it the 
well-known coefficient of the anomaly. 
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We may now verify this by evaluating the two-gluon matrix element of the 
renormalized anomaly equation ( |3.14| ) using the results obtained above. There are 
no contact terms. The reduced matrix element of d^[j^]Qg — 2m[j5]os was found 
to be the same in the two kinematical settings studied through two loops, 

But this is exactly equal to /3^'* ^[ff]ms limit p' ^p, so the validity of the 

renormalized anomaly equation requires that there be no oversubtraction for FF, 
i.e., that [FF]ms and [FF]Qg be the same object. 

The verification of this consistency requirement, which previous treatments 
P,pil| have omitted altogether, occupies the remainder of this section. One way of 
accomplishing it would be to directly compute the term of order g"^ in fs that gives 
the finite renormalization relating [FF]ms to [FF]Qg at second order. I will forego 
this option for a less direct and therefore maybe more interesting check. [FF]Qg 
is characterized among the finite multiplicative renormalizations of [FF]ms by the 
property that g^[FF]Qg has no anomalous diagonal scaling. Ergo, [FF]ms is equal 
to [FF]Qg if and only if its diagonal anomalous dimension is exactly (3{g'^)/g'^. 

For this to be the case the minimal subtraction constant Z^p must be equal 
to Z^^ through two loops; we already know that to be true at one loop. In view of 
( ^77| ) this amounts to i'pp being finite if its parameters are expanded to order 

The thirty-six two-loop diagrams contributing to t'pp ^^^^ are shown in Figure 
7. Diagrams in which both external gluons couple to the same four-point vertex 
may be discarded immediately by symmetry considerations and have been omitted; 
no orientations are indicated on fermion and ghost lines and it is understood that 
both choices are to be considered. Table 4 lists the pole parts of these diagrams 
in the Feynman gauge (a = 1) as computed for p' ^p; for the diagrams involving 
fermions, the limits p^/m^— >0 and m^/p^— >0 are distinguished if necessary. Each 
contribution of the form ^ is to be supplemented by a —^Lp term of equal strength, 
left out of the table for the sake of brevity. 
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The contributions of order carry three distinct group theoretical structures, 
namely c\, CAnnTu, and CfiUpTfi. The cancellation of infinities we are after should 
obtain for each structure separately. 

Consider first the c\ terms. Retaining only those terms is tantamount to 
assuming that the fermions are absent from the theory, which reduces the set of 
diagrams to numbers I through XXIX. The sum of these diagrams as per Table 4 is 

The cancellation of ^ poles in this sum, with the accompanying cancellation of 
terms of the form ^Lp, reflects the absence of divergences in the one- loop matrix 
element [B^. The counterterms follow from expanding Qq and the bare gauge pa- 



rameter ckq in ^ FF barelone loop in reuormallzed quantities; if the beta-function is 
written as big"^ + . . . and the renormalization group coefficient of the gauge param- 
eter, H-j^cx, as dig"^ + . . . , that amounts to applying 

e dg'^ e da 



to the expression in ([4.22|) . With the well-known values bi = — llc^/247r^ and 
di = 5c^a/247r^ that leads in the Feynman gauge to a counterterm contribution 



exactly opposite to (|4.28| ) . This shows that the g^ part of Tpp is finite and 



hence to this order in pure QCD Zpp = Z^^, as desired. 

Upon re-inclusion of the fermions there are seven more two-loop diagrams, each 
involving a single fermion loop, shown as numbers XXX through XXXVI in Figure 
7. The Feynman gauge pole parts of these diagrams, again computed at p' =p, are 
sensitive to the value of /w? . Adding up the relevant entries of Table 4 produces 



I Jlc^nnTa \^A^^!^ ifmVp^-0; 

Concomitantly the values of hi and di receive the well-known fermionic contribu- 
tions 5hi = ubTr/Q-k'^ and 5di = —riRTfia/GTv'^. It is easy to check that the pole 
terms with group structure caTr cancel between the new diagrams and the {a = l) 
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counterterm modifications. Thus it only remains to account for the cunpTji part 

Inspection of Table 4 shows that the net CRnnTji contribution arises from di- 
agram XXXI. The upper part of that diagram is identical to the one-loop mixing 
diagram of FF and a^j^, Figure 3. In fact it is not FF but FF -f ^^d^j^ that 
is multiplicatively related to [FF]mSj where the pole is that of the mixing diagram; 
this formula accounts for (|2.7|) , with In ^ replaced by j . The renormalization for- 
mula for FF, ( |4.7D , should be modified accordingly. In view of the one- loop values 
of Tg^ji^ij)' = p) in the two limiting regimes of /m? , (4.5-6), this effect provides a 
correction that precisely cancels the Ci^ri/jTR-dependent term in (|4.30|) . The pecu- 
liar dependence of ( 4.30 ) on the presence or absence of a fermion mass is thus seen 
to be inherited from the analogous feature of the divergence of the axial current. 

This concludes the consistency check (at least as far as the Feynman gauge is 
concerned): through two loops the minimally subtracted operator [FF]ms scales 
as the inverse of the minimally subtracted coupling constant, without the need for 
a finite renormalization. 



5. Comments 

In the previous section the two-gluon matrix elements of a number of suitably 
normalized operators were computed through two loops in order to verify the valid- 
ity of the renormalized anomaly equation. I will now consider those matrix elements 
in their own right. 

The first issue that should be addressed is the nature of the various kinematical 
limits. Most of the explicit results listed above have applied to the reduced matrix 
elements of composite operators at total momentum p' — p equal to zero, leaving 
the momentum squared of each external gauge field, p^, and the mass m of the 
fermions as independent parameters. Existing physical applications typically involve 
the reduced matrix elements in the limits in which —p'^ goes either to zero or to 
infinity with respect to m^. It is well-known that at one loop the reduced matrix 
element of the current divergence vanishes as — p^/m^ — s>0, whereas the pseudoscalar 
density mjs drops out in the opposite limit, see (4.5-6). Evidently neither of these 
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observations depends on the precise normalization of the objects they apply to. 
Both statements remain true at two loops, the first, exhibited in (|4.10| ), being 
obviously the non-trivial one. Such behavior has been observed in the formalism of 
|T0[] as well. The noncommutativity of the infrared limits has also been discussed 



at the two-loop level in a dispersive approach in |^^; however, the contention of 



that paper that different limits correspond to different currents is at odds with the 
point of view taken here. 

For the abelian theory the vanishing to all orders of the reduced matrix element 
of d^j^ in the limit — p^/m^ — >0 is an aspect of the original Adler-Bardeen analysis 
0; the proof uses gauge invariance and analyticity. The background field method 
transfers at least the former tool to the nonabelian case. A putative generalized 
proof would start with the observation that if symmetry considerations are taken 
into account, the IPI function of the (renormalized) current and two gauge fields 
has the form [^, at p'"^ =p^, 

(5.1) 

+Ds{p^) {p^e^^pp'^p'^ + py.c.f)P"p'^)\ Sab 

and the (renormalized) reduced matrix element of dxj^ is —2iDi. Gauge invariance 
of the background field effective action requires ( p.l|) to vanish when contracted with 

QY p'^ ^ implying DiijP') = — p^(L>2(p'^) + -D3(p^)). If -D2 and are sufficiently 
well-behaved at = it follows that Diijp' = 0) = 0; but the premise is not at 
all obvious in the nonabelian case, whereas the L>'s are known to be analytic near 

= in the abelian theory with massive fermions. 
The interpolation between the two extremes —p^jw? — and m? /{—p^) -^0 
involves m explicitly and is therefore highly scheme-dependent. At low orders, 
however, relatively simple results may emerge; for instance, for the two- loop (1 — a)^ 
terms of ^df^[j'^]Ms -^"i[j5]Ms ^'^^ interpolation turns out to exactly parallel that 
of the complete one-loop functions, as in ( |4^ ) and ( f4.3|) . 

Granted the distribution of nonzero results between the divergence of the cur- 
rent and the pseudoscalar density, what are those results taken by themselves? 
The simplest case is the massive abelian one at p = p', ^ 0, treated by Adler 
and Bardeen 0|. In that case the anomaly equation relates the matrix element 
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of the canonical pseudoscalar density m[j5]c to that of g'^[FF]c; the latter is no 
longer sensitive to mixing because the matrix element of d^j^ vanishes under the 
circumstances. (In this sense what is often referred to in the literature as the Adler- 
Bardeen current is merely an accessory to the Adler-Bardeen theorem.) It can 
be shown by inspection of diagrams that the complete reduced matrix element of 
g^[FF]c, hence also that of mfisjc, is given by its lowest order value, a constant 
times (7^. Each of these two reduced matrix elements satisfies a renormalization 
group equation 

+ ^ir2 + ^ + /^^TT - 27A Toip' =P,P^ = 0;m^0) = (5.2) 

where the anomalous dimension of the photon is essentially identical to the beta- 
function, 7yi = /?/2(7^; the exact results for the Toip' =P',p'^ = 0;my^O) thus satisfy 
the renormalization group equations ( |5.2|) in the minimal nontrivial way |jTl|], with 
IJi-^Vo = 0. The explicit two-loop values, (^4.19|) and (|4.22|) at c^ = 0, are obviously 
compatible with the general result. 

Still in the abelian theory but with vanishing fermion mass, it becomes im- 
possible to disentangle the divergence of the axial current and the antisymmetric 
photon tensor by the infrared behavior of their reduced matrix elements, because 
of mixing and the attendant normalization ambiguity. Perhaps the cleanest matrix 
element in this kinematical configuration is that of the divergence of the canonical 
axial current, dfj\j'^]c-, which satisfies a renormalization equation of the form ( |5.2|) . 
However, the vanishing anomalous dimension of [j^]c now merely conceals the un- 
derlying non-trivial renormalization effects and the solution to ( J5.2D is no longer 
exhausted by a one- loop term. Explicitly, from ( |2.12| ), ( |2.11| ), ( p.7| ), and ( [4.17| ) 
(with Cyi = 0), and using the lowest order value of /3, 

^d,[j^]ciP =P;m = 0)\abeHan = ^4 +■••• (^'S) 

The renormalization group calls for a logarithmic term at order g^ to match the g'^ 
contribution; this term is supplied by three-loop diagrams as in Figure 8, computed 

2 

and discussed in [|2^ . Of course higher powers of In will emerge from the ellipsis 
if the calculation is pursued to higher orders. 
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Turning now to the nonabelian case, consider only the case of massive fermions 
and p' = p,p'^ —>-0, obviously the simpler. In ordinary perturbation theory the anal- 
ogy with the abelian situation is obscured from the onset, because the anomalous 
dimension of the gluon, 7^, is then no longer simply related to the beta-function; 
this spoils the chances of a purely lowest order reduced matrix element, as this 
could no longer saturate the renormalization group equation. In the background 
field method, 'Ja is again equal to ^ p-. As a result there can be no logarithms in the 
reduced matrix element of mfjsjc' or g^[FF]c at the next-to-leading order, a fact 
borne out by the explicit results, (|4.19| ) and ( [4.22|) (in contrast to the corresponding 
quantities in [jl^, obtained in ordinary perturbation theory). Yet, as these equa- 
tions show, even in the background field formalism r^[jg]^ (p' =p, — > 0; m 7^ 0) and 
g2[pp^^{p' = Pt p'^ 0', m 0) have corrections to their leading orders, in contrast 
to the abelian case. 

A salient feature of these corrections is their dependence on the gauge pa- 



rameter a. That an expression for general p and p' such as ([4.20D exhibits such 
dependence is not surprising, because there is no reason to expect gauge invariant 
results as long as the external gluons are off-shell. Of course infrared problems make 
it impossible to take p^ —*0 or p'^ -^0 in ( [4.20|) . If p' is first taken to be equal to p 
the one-loop result becomes independent of and p'^, see (|4.22| ), and in particular 



finite as p^ —>■ 0, but this is deceptive; the vanishing of the anomalous dimension 

of g'^[FF]c is accommodated in r^2[F^]^(p' = p) by the occurrence of powers of 

_ 2 

In starting at the two-loop level. 

The noncommutativity of these limits is reminiscent of what happens with the 
divergence of the axial current at one loop. The analogy may be reinforced by 
introduction of a small finite gluon mass M. This was probably first considered by 
Chanowitz [^^, who employed spontaneous symmetry breaking to procure a mass 
for the gauge fields. As long as gauge fields contribute through no more than one 
loop the mass may also be introduced straightforwardly, via a massive propagator 



(5.4) 



.k'^-M^ {k^ - M^){k'^ -aM^) 
Recalculation of the diagrams in Figure 2 with this propagator leads to a general 
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ization of (|4.20| ). In the Feynman gauge, 

i'lFFuJone loop = -^^^(l + (2^^ " ^ip - p'f) hoiP: p' ^ M)^ (5.5) 

This should be compared with the expression for the one-loop reduced matrix el- 
ement of d^j^, ([4.2|) . At p' =p, i.e. as the composite operator is integrated over 
spacetime, the two results are completely analogous. In particular, as long as M is 
finite the one-loop correction to i'^pFjMs disappears at p' =p. This is actually also 
true for gauges a 7^ 1, though the expression for general p and p' is not quite as sim- 
ple as ( |5.5| ) above. The impossibility of a direct generalization of the Adler-Bardeen 
results to the nonabelian case appears therefore attendant upon the masslessness of 
the gauge fields. 

We may reconsider similarly the calculation of iL[j^jj^jg{p=p' ,p'^ ^0; rriy^O) by 
including a small gluon mass, taking the limit ^ — * 0. Not surprisingly, the entries 
in Table 3 that betray their sensitivity to the infrared behavior of the gauge fields 
by the presence of ln(— p^) have to be changed to reflect the new circumstances. 
Diagrams I through VI, VIII and XI are not affected, but VII now yields 

(-12-4a)ln— -4ctlna + 2 + 6a (5.6) 
V / 

and IX+X 

/4a \ 

2ca[ 4aL^ + (6 + 2a) In — + 2a In a - 2 - 4a . (5.7) 

V e / 

Compared with table 3 as it stands the extra contributions add up to — 16c^[l -|- 
|(1 — a) — -1(1 — a)^], which is precisely what is needed to eliminate the terms 
proportional to 0^(7*. 

The existence of a finite one-loop contribution to the reduced matrix element 
of g'^[FF]c aX p' = p appears to be at odds with a semiclassical interpretation 
of ■Q^g'^[FF]c as the Pontryagin density, the integral of which is the instanton 
number. It is tempting to try to remove the extra factor by a finite renormalization, 
but that would interfere with the vanishing of the anomalous dimension of [FF] c 
and leave explicit logarithms untouched. However, the sensitivity of the one-loop 
reduced matrix element to the infrared arrangement of the theory suggests that 
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the problem is an infrared artifact of perturbation theory, and that its vanishing at 
infinitesimal gluon mass mimicks behavior induced by a nonperturbative physical 
infrared cutoff. (This appears to be also the point of view advocated in [|15|; the 
alternative point of view that all or part of the one- loop matrix element of g^FF 
is a signature of a real effect involving the vacuum angle 9qcd [HI is developed in 



6. Conclusion 

Since its identification more than twenty years ago, the axial anomaly has never 
been far from the spotlight of theoretical physics. Regarding most of its aspects 
it could therefore be plausibly argued that no relevant insight remains to be put 
forward. Yet the continuing fiow of papers related to the interplay of the anom- 
aly and the process of renormalization bears witness that the concrete organization 
of the currently available intuition in specific physical contexts is still not with- 
out challenges. The rationale of the present paper was the consideration that for 
nonabelian gauge theory and the practically important formalism of dimensional 
regularization, a truly comprehensive discussion was lacking in the literature. 

A review was given of the renormalization group arguments that lead to a 
convenient normalization of the singlet axial current and the operators related to 
its divergence. These considerations were matched to a derivation of the anomaly 
equation using dimensional regularization with what is probably the most straight- 
forward treatment of pseudotensorial objects, that according to 't Hooft and Velt- 
man. This settled the issue in principle. To illustrate how this works in practice, 
and verify that the abstract argument was sound, I computed the two-gluon ma- 
trix element of each term in the anomaly equation through two loops, including 
features neglected in previous treatments. Such calculations supply, apart from the 
desired reassurance, matrix elements that may be of independent interest, if not of 
immediate physical significance. 
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Table Captions 



Table 1. Two- loop contributions to the unrenormalizcd reduced matrix element of f^^jg 
at p' =p and m = 0, numbered as in Fig. 5. A common factor g^URTu/ {A.t^Y is 
omitted. 

Table 2. Two-loop contributions to the unrenormalized reduced matrix element of f^^jg 
for p' = p, p'^ and m ^ 0, numbered as in Fig. 5. A common factor 

g^niiTji/{4n)^ is omitted. 

Table 3. Two-loop contributions to the unrenormalized reduced matrix element of js 
for p' = p, ^ and m ^ 0, numbered as in Fig. 5. A common factor 

g^nBTn/ {ATx)^m is omitted. 

Table 4. Pole parts of the two-loop Feynman gauge contributions to the unrenormalized 
reduced matrix element of FF for p' = p, numbered as in Fig. 7. A common 
factor g^/{4Tv)'^ is omitted, all terms of the form ^Lp. In XXXI, (*) 

marks the result for m — 0, p^y^O, (**) that of rriy^O, p^^O. 



Figure Captions 

Fig. 1. One- loop diagram relevant for the renormalization of and j^. 

Fig. 2. One loop diagrams contributing to the two-gluon IPI function of FF. 

Fig. 3. One- loop diagram leading to mixing of d/^j^ with FF. 

Fig. 4. One- loop diagram for the two-gluon IPX functions of d^j^, J5, or Oan- 

Fig. 5. Two-loop diagrams contributing to the two-gluon IPX functions of dfj^j^ and 

Fig. 6. One- loop diagrams contributing to the ip-ip IPX function of Oan- 

Fig. 7. Two-loop diagrams contributing to the two-gluon IPX function of FF. For each 
ghost or quark line both orientations should be taken into account. 

Fig. 8. A three-loop contribution to the two-gluon IPX function of d^j^. 
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TABLE 1 



(l-a)° 



I + II 2cr(^^ -12Lp + 2l) 2cr(^-^ + 12Lp-2l) 



IV+V 2(cr - |ca) (-^ + 8Lp - is) 2(cr - ^ca) - SL^ + m) 



16 _\ /16 







III cr(^-8Lp + 14) cr(-^ + 8Lp-14) 







VI cr(^-— + 16Lp-52) cr(^— - 16Lp + 28) 

VII CA^-y + 32Lp - 104) ca(^-8Lp + 10) 

VIII 
IX+X 2CA - 4Lp + y + 8C(3)) 2cA - 2Lp - ^ - 4C(3)) 4ca 

XI CA(y-16Lp + ^-16C(3)) ca(-^+4Lp-11 + 8C(3)) 

total -32cr-32ca -24ca 4ca 



TABLE 2 



(1 - af 



(l-af 



I + II 
III 



VI 

VII 

VIII 

IX+X 

XI 



/4 ii 
2cR - -4L^ - — , 
\e 3 / 

/ 24 10\ 
Cr(^-— + 24L^ + yj 

IV+V 2(cr - |ca) - 8L^ + 2) 



/ 32 \ 
ca(-— + 32L^ + 8j 



2ca(— -12L„-l) 



ca(^Y ~ ^^^^ 



2cr(^-4L^-3) 
Cr(^ -8L„-2) 
2(cR-|cA)(-- + 8L^ + 2j 



e 

4cr 



Ca(^ -8L^-2) 



2ca(--+6L„ + -) 
ca(-^ +4L^ - 1) 















total 















TABLE 3 



(l-a)° 

I + 11 2cr(^ + 16L^- y) 2cr(^-8L^-2) 

III Cr(-^ + 16L^- y) Cr(^-8L„-2) 

IV+V 2(cr-|ca)(^-8L^ + 2) 2(cr-|ca)(-^ + 8L^ + 2) 

VI CR(y-32L^ + 20) cr(-^ + 8L^ + 2) 

VII ca(-16Lp + 16L„ + 56) ca (4Lp - 4L^ - 

VIII 
IX+X 2cA + 8Lp - 12Lm - 22) 2ca (- - - 2Lp + 6Lm + lo) -2ca 
XI 6cA 

total 8cr + 16c A 12ca — 2ca 



TABLE 4 



I 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 

X 

XI 

XII 

XIII 

XIV 

XV 

XVI 

XVII 

XVIII 



148 541- 



320 



1 4 



1104\ 2 



/16 2\ 2 

/12 35\ 2 



(- 




608 


5960- 


3e2 


9e , 


304 


2716 


3e2 


9e . 


32 


344 \ 


~3? 


9e / 


16 


196 \ 




9e / 




18 2 



/ 62 
V 3? 



/144 396 \ 







/108 243 \ 2 

1^ + — jci 



V e 



192 


+ 


576 ^ 








)ci 


£2 








288 




792 ^ 




e2 


+ 




)ci 



XIX 

XX 

XXI 

XXII 

XXIII 

XXIV 

XXV 

XXVI 

XXVII 

XXVIII 

XXIX 

XXX 

XXXI 

XXXII 

XXXIII 

XXXIV 

XXXV 

XXXVI 



48 2 

192 624 
— + 

192 



4 



480 \ 



108 



255 \ 2 
144 956 \ 2 



+ 



3e 



■)' 



/16 140\ 2 



192 
36 



384 



ci 



192 432 



-A 



^CAnRTn 



128 ^ , , 64 ^ , , 

CrHrTr (*), CrHrTr (**) 



64 



e 

512 



(cr - |cA)nRTR 



,3e2 ^ 


9e . 


256 


2176- 


,3.^ + 


9e , 


/64 




t + 






e 


128 


800 


£2 


3e 



R 



jcAHRT 
^CaHrTr 

)cAnRTR 



